D'yakonov-Perel' spin relaxation in InSb/AlInSb quantum wells 
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We investigate theoretically the D'yakonov-Perel' spin relaxation time by solving the eight-band 
Kane model and Poisson equation self-consistently. Our results show distinct behavior with the 
single-band model due to the anomalous spin-orbit interactions in narrow band-gap semiconductors, 
and agree well with the experiment values reported in recent experiment (K. L. Litvinenko, et al., 
New J. Phys. 8, 49 (2006)). We find a strong resonant enhancement of the spin relaxation time 
appears for spin align along [110] at a certain electron density at 4 K. This resonant peak is smeared 
out with increasing the temperature. 

PACS numbers: 72.25.Rb, 71.28.+d, 71.70.Ej. 



i 

C 

o 
o 



> 
in 
in 

00 

o 



OS 

o 



X 



Spin relaxation time (SRT) is very important for the 
coherent manipulation of electron spin and applications 
in spintronics devices. There are four different spin re- 
laxation mechanisms^, i.e., the D'yakonov-Perel' (DP) 2 , 
Elliott- Yafet (EY) 3 , Bir-Aronov-Pikus 4 and hyperfine 
interaction^ mechanism. Among them, the DP mech- 
anism is found to be the dominating spin relaxation 
mechanism in zincblende semiconductor structures over 
a wide range of temperature^. According to the DP 
theory, the electrons lose their initial spin orientation 
due to a momentum-dependent effective magnetic field 
that changes its orientation frequently which is caused by 
random impurity scattering. Therefore the momentum- 
dependent effective magnetic field is the key factor to de- 
termine the spin relaxation time, and it is induced by two 
types of spin-orbit interactions (SOIs) in structures with- 
out inversion symmetry, i.e., the Rashba SOI (RSOI) 
arising from structure inversion asymmetry and the Dres- 
selhaus SOI (DSOI) caused by bulk inversion asymme- 
try. In conventional semiconductor quantum structures, 
the DP theory based on the single-band model with lin- 
ear momentum-dependent SOIs have been demonstrated 
to agree well with the experiments, e.g., GaAs/AlGaAs 
quantum well (QW) 7 -, and InGaAs/InP QW— . Recently, 
the spin relaxation time in narrow band-gap semiconduc- 
tor InSb/AlInSb QW also attracted much interest ^ 10 ' 11 , 
because of its unusual properties for spintronics devices, 
e.g., small electron effective mass, strong spin-orbit cou- 
pling, large effective Lande g factor. However, the SOIs 
in narrow band-gap semiconductor is quite different from 
the single-band model with linear momentum-dependent 
SOIs. For example, the Rashba spin-splitting exhibits a 
nonlinear behavior while the kinetic energy of the elec- 
tron is comparable to the band-gap^. Therefore a de- 
tailed theoretical investigation beyond the single band 
model for the DP spin relaxation time is necessary in 
the case of narrow band-gap InSb/AlInSb QW for the 
potential spintronics device and basic physics. 



electron density and the thickness of the InSb/AlInSb 
QW. The effective magnetic field and the spin relaxation 
time is calculated based on the self-consistent solution of 
the eight-band Kane Hamiltonian and the Poisson equa- 
tion. We find that the effective magnetic field obtained 
from the eight-band model deviates strongly from that 
obtained by the single-band model with the momentum- 
linear SOI. We show that the eight-band model results 
are in good agreement with experiment values without 
introducing any fitting parameter. We find a strong 
anisotropic SRT: a strong resonant peak of SRT for spins 
aligned along the [110] direction can be seen by tuning 
the electron density since electron at T = 4 K. But this 
peak is gradually smeared out with increasing tempera- 
ture. Our results could be helpful to observe new physical 
phenomenon, e.g., the intrinsic spin Hall effect and per- 
sistent spin helixi 3 - in such narrow band-gap InSb/AlInSb 
QWs. 

We consider an asymmetric n-doped InSb/AlInSb QW 
grown along the [001] crystallographic direction (see Fig. 
Q] (a)). The n-doping layer is assumed to be located 20 
nm on the left-side of the InSb well and with an expo- 
nentially decaying profile. We extend the previous theory 
(see Refs. [TH andllot) to the framework of the eight-band 
model by changing all operators in two-band model to 
the eight-band mode l 17 i 18 , the DP spin relaxation time 
r a (a = +, — , z, representing the spin relaxation time for 
the spin of the injected electrons oriented along [110], 
[110], [001], respectively) can be written as 
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In this work, we investigate theoretically the DP spin 
relaxation time and its dependencies on the temperature, 
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Here, E^(k) = E s (k) — E s (0) is the kinetic energy 
of an electron in the s-th subband, AF Si+ (Epk) = 
F s%+ (EF.k) — F s ,_ (Epk) is the Fermi distribution dif- 
ference between the spin-up and spin-down subband, Tt r 



is the transport relaxation time and r s>a (fe) (a 
is the spin relaxation rates 



r Si+ (fc) = A XX + Ayy - A X y ~ A 
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FIG. 1: (a) Fhe calculated band alignment and electron (hole) 
probability distribution. The shading area describes the dop- 
ing profile. The sheet carrier density is n e = 5.2 x 10 11 cm -2 , 
(b) The self-consistently calculated energy dispersions of elec- 
tron, (c) The vectors of effective magnetic field fli(fe) on 
k x -ky plane (d) The magnitude of the effective magnetic field 
fii(fc) for k along [100] (black solid line) [110] (blue dashed 
line) and [110] (red dash-dotted line) of an asymmetrically 
n-doped InSb/Alo.isIno.ssSb QW at T = 4 K. 
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^s.i(k) (i = x, y) is the components of the in-plane effec- 
tive magnetic field of the s-th subband. f2 s ^(fe) can be 
obtained by ascribing the spin-splitting induced by space 
inversion asymmetry to the Zeeman splitting caused by 
the effective magnetic field. Therefore, by using the 
eight-band Zeeman termi£ H z = /i^B • S, J2 s> i(fe) can 
be written as 



O s ,i(fe) = fisBi 



S s 4k)AE s (k) 



S s ,x(k) 2 + S s ,y(k) 



2 - 



(10) 



S s .i(k) = (^ s ,+ (fe)|E i |V s , + (fe)> - 



-(*)>• 
(11) 

£,(i = x,y) is the components of the eight-band effec- 
tive spin matrices which can be found in Ref. [TH and 
[H. AE s (k) ee E Si+ (k) - E s _(fe) is the spin-splitting of 
the s-th subband. The eigen-energy E St ±(k) and eigen- 
states |^ Si -|_(fe)) can be numerically obtained by solving 
the eight-band Kane Hamiltonian and the Poisson equa- 
tion self-consistently 18 . Through this approach, the non- 
parabolic effect and the anomalous behavior of SOls in 



narrow band-gap semiconductors can be taken into ac- 
count. In Eq. (JTJ)- (JSJ , v is a constant characterizing the 
relation of momentum scattering time on the electron 
kinetic energy (r p (fc) oc [E r s (k)Y). For acoustic phonon 
and screened ionized impurities scattering (type I), v = 0, 
for polar optical phonon scattering (type 11), v = 1, for 
weekly screened ionized impurities (type III), v = 2i£. 

In Figs. HJa) and (b) we show the calculated potential 
profile, the electron probability and the energy dispersion 
of a 10 nm InSb/Aln.isInn.ssSb QW at T = 4 K including 
the effect of built-in electric field caused by the charge re- 
distribution. All the Kane parameters of the materials 
used in our calculation are taken from Ref. [20L and the 
ratio of the conduction band offset and the conduction 
band offset is taken as 62%:38% 2 i. The bulk inversion 
asymmetry of zincblende crystal is introduced by the B 
parameter in the eight-band Kane Hamiltonian 22 , which 
is taken to be B — 31.4 eV ■ A 11 . Besides, we should 
notice that the temperature dependence of the bandgap 
(Varshini relation)^ is more pronounced in a narrow 
bandgap semiconductor than that in a wide bandgap 
semiconductor, for instance, the bulk bandgap of InSb, 
of which the bulk band gap is 0.235 eV at 4 K, and 0.174 
cV at 300 K, i.e., up to a 26% variation of the bandgap 
with increasing temperature. Because the SOI and spin- 
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FIG. 2: Calculated spin relaxation time r z versus tem- 
perature compare to the experimental results reported in 
Ref. |9|. (a) Sample 1833 (asymmetrically n-doped 20 nm 
InSb/Al .i 5 In .85Sb QW) (b) Sample 1831F (uniformly n- 
doped 20 nm InSb/Al .i 5 In .85Sb QW). The black, blue, red 
lines represent the eight-band numerical results for type I 
(y — 0), type II (y = 1), and type III (y = 2) momentum 
scattering mechanism. The dashed lines with the same colors 
represent the results from the single-band model with linear 
SOI. 



splitting is intimately related to the conduction-valence 
band coupling, the decreasing of band gap could lead to 
an enhancement of the SOI and results in the increas- 
ing of electron subband spin-splitting and the effective 
magnetic field (about 11%). 

Fig- HJc) shows the effective magnetic field as a func- 
tion of the in-plane momentum. The effective magnetic 
field of a (OOl)-grown InSb/AlInSb QW always lie in the 
QW plane. Due to the interplay of RSOI and DSOI, the 
effective magnetic field exhibits a Civ symmetry. Notice 
that the effective magnetic field is enhanced or weak- 
ened when k along the [110] or [110] directions. Fig. 
[]Jd) shows the self consistent eight-band modeling for 
the magnitude of the effective magnetic field for k along 
[100], [110], [110] (see the solid curves). Along [110] crys- 
tallographic direction, the effective magnetic field point- 
ing along [110] vanishes at a certain Fermi wavevector k, 
which makes the spin lifetime t_ become very long. In 
Fig. Q] (d) we compare these results with those from the 
single-band model with fe-linear SOI. One can see clearly 
that the deviation of effective magnetic field is very large 
(up to 42%) at k = 0.2 nm" 1 ([100]). This deviation 
comes from the weakening of the interband coupling as 
the electron kinetic energy becomes comparable to the 
bandgap. Therefore the single-band model may not be 
good enough to describe the strong SOIs in such narrow 
bandgap QWs. 

In Fig. [5] we compare the numerical results of SRT t 2 
with the experimental measurement reported in Ref. [9|. 
The transport momentum relaxation time r tr used in our 
calculation are obtained from the measured Hall mobil- 
ity (by T tr = m*fJ,Hall/ (erHaii)^)- The electron density 
is assumed to increase linearly from 3.6 x 10 11 cm~ 2 to 
5.3 x 10 11 cm" 2 for sample 1833 and from 5.7 x 10 11 




FIG. 3: Spin relaxation times t z , r+ and t_ as a func- 
tion of electron density n e in a 10 nm asymmetric n-doped 
InSb/Alo.i5lno.s5Sb QW for different temperature. The 
dashed lines with the same colors represent the results of 
single-band model with linear SOI. 



cm~ 2 to 7.3 x 10 11 cm -2 for sample 1831F when temper- 
ature increase from 77 K to 300 K. Considering differ- 
ent momentum relaxation mechanisms, our eight-band 
numerical results agree quite well with the experiment 
values without having to introduce any fitting param- 
eter. From panel (a), we can see the weekly screened 
impurity scattering (y = 2) and polar phonon scatter- 
ing {y = 1) dominate at T < 150 K and the ionized 
impurity scattering dominates at T > 150 K. Noticed 
that in heavily doped semiconductor samples, the dom- 
inant momentum scattering mechanisms varied through 
neutral (weekly screened) impurities scattering, acous- 
tic and polar phonon scattering, and ionized impurity 
scattering with increasing temperature^, therefore our 
results are reasonable and also consistent with the pre- 
vious wor k 16 i 24 . One can see that, due to the overesti- 
mate of SOI strength by the linear SOI model, the single- 
band model will underestimate the SRT compared to the 
eight-band model and doesn't agree with the measured 
SRT. For the uniformly doped sample 1831F, our calcu- 
lated DP SRT is larger than the measured value. The 
discrepancy between the calculated DP SRT from the 
eight-band model with measured value is reasonable be- 
cause the SRT induced by the EY mechanism could be 
comparable to the DP SRT 9 . The calculated SRT from 
the single-band model in this symmetric doped sample 
is similar with that from the eight-band model. This is 
because the cubic DSOI term (the RSOI is absent due to 
symmetric doping) in this sample may play a dominant 
role. 

In Fig. [3] we calculate the spin relaxation times r z , 
t+ and t_ as a function of electron density n e in a 10 
nm n-doped InSb/Alo.isIno.ssSb QW for different tem- 
peratures. As shown in the figure, the SRTs r 2 , t + de- 
crease with increasing electron density due to the en- 
hancement of SOIs with increasing the Fermi wavevec- 
tor. When temperature increases, the r z and r_ are 
suppressed strongly, but r+ is not very sensitive to tern- 
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FIG. 4: Spin relaxation times t z , r+ and r_ as a function 
of carrier density in asymmetric n-doped InSb/Alo.isIno.ssSb 
QWs. The black solid, blue dashed, red dashdot and green 
dotted lines represent the numerical results of Li„sb = 
10, 15, 20 and 25 nm. The dashed lines with the same col- 
ors represent the results of the single-band model with linear 
SOI. 



perature. The giant spin relaxation anisotrop} ! 11 ! 15 ] 16 is 
also demonstrated in InSb/AlInSb QW: the SRT of the 
[110]-oriented spins shows a resonant peak at a certain 
electron density at low temperature. However, the reso- 
nant peak is very sensitive to temperature. As temper- 
ature increases from 4K to 300K, the peak is gradually 
smeared out, which is due to the bluring of the Fermi 
surface with increasing temperature. The dashed lines 
show the results of the single-band model. Consistent 
with the former discussion, the results of the single-band 



model is smaller than that of the eight-band model but 
at small n e . 

In Fig. [4] we exhibit the SRTs r z , r+ and r_ as 
a function of electron density in a 10 nm n-doped 
InSb/Alo.isIno.ssSb QW for different thicknesses of InSb 
well. The SRTs increase with increasing thickness of the 
InSb well when the electron density n e is small, but has 
an opposite trend when the electron density n e are larger 
than a certain value. The resonant peak value of r_ be- 
comes larger and the corresponding n e becomes smaller 
with increasing thickness of the QW. The results of the 
single-band model is smaller than that of the eight-band 
model except at small n e . 

In summary, we investigated theoretically the SRT in 
InSb/AlInSb QW beyond the single-band model. Our 
results are obtained within the eight-band model and 
agree very well with the measured SRTs, while the 
SRT obtained from the single-band model with linear 
momentum-dependent SOIs deviates strongly from that 
of the eight-band model due to the strong interband- 
coupling in narrow bandgap QWs. We also demonstrate 
that the SRT along [llO] -direction shows a resonant peak 
at a certain electron density, i.e., very long SRT. The res- 
onant peak will be smeared out with increasing temper- 
ature. 
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